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ABSTRACT 



There is a need for this country to bo able to assure itself that 
an international ban on nuclear testing can bo made effective* This 
thesis studies one small portion of that problem - to wit, the procedure 
for classifying an underwater detonation which has vented and which has 
been detected by a hydroacoustic network* 

A MEASURE OF EFFECTIVENESS is developed, and after a st\idy of some 
parameters such as height of search, search spacing, search width and 
current movements optimizing search plans are deduced* 

One of the major difficulties in such a quantitative analysis is 
the non-doterministic nature of many of the parameters* The ideal 
search spacing was found to depend directly on the surface radioactive 
pool diameter* In order to develop plans this essential factor must be 
known* thifortimately this depends on many variables over which the 
searcher has no control or knowledge* An averaging procedure is adopted 
and a factor established to compeure spacings and determine which meets 
the criterion dictated by the MEASURE CF EFFECTIVENESS* 
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1 . Introduction. 

In October 1960 a technical sympoaium was held by the Advanced 
Research Projects Agency (ARPA) acting in behalf of the Department of 
Defense for the purpose of informing those in attendance of the scope, 
status and planning underway for Project VELA. Project VELA consists of 
research and development that is aimed at the improvement of our capa- 
bility of detecting and identifying nuclear detonations that are con- 
tained below the siirface of the earth, or underwater or that occur in 
space above the atmosphere. 

It is well recognized that it is p>ossible for a clover and deter- 
mined nation to conduct subsurface and high altitude nuclear weapons 
tests without a significant risk of being detected. Detection of 
nucleair detonations on the surface or within the atmosphere is consid- 
ered well in hand and not included in VELA. Ever since the Conference 
of Escparts held in 1958 in Geneva, international group® have been 
attempting to establish methods and networks which would assure the 
participating nations that agreements to bfui nuclear testing could be 
effectively enforced and thereby be the first step in solving one of the 
most important paroblems of our time - ARMS CONTROL. 

The hope of Project VELA is to improve our detection and classifi- 
cation capability by at least one order of magnitude, and parhapa by as 
much as three, over a three year span. Such a quantum jump is considered 
feasible if sufficient independent scientific interest can be aroused 
while assigning a high priority to the investigation in as many of the 
government laboratories which have some measure to contribute. 

The Project was subdivided into three major categories: 
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VELA HOTEL - Detection of very high altitude nuclear detonations by 
means of detectors mounted in satellites; 

VELA SIEPIRA - Development of techniques and equipment to detect and 
identify high altitude nuclear detonations by ground based instrumenta- 
tion; 

VELA UNIFORM - Improve, markedly, equipment and techniques for de- 
tecting and identifying subsurface nuclear explosions from a distance as 
well as to develop rapid, economical techniques for conducting on-site 
inspections to determine whether a suspicious event was actually a nuclear 
explosion or not* 

Project VELA UNIFORM includes a study of underwater detonations as 
part of its ”subsurface” objective* The Naval Electronics Laboratory 
(NEL) at San Diego was assigned the task of investigating the detection 
of underwater explosions and the Naval Radiological Defense Laboratory 
(NRDL) was contracted to devise a searching procedure for obtaining sam- 
ples to classify explosions detected by a network devised by NEL* 

In a CONFIDENTIAL report /6/ NEL outlined a distribution of hydro- 
acoustic stations which would be capable of detecting underwater blasts 
in most areas of the world* However, this network has practically no 
ability to classify explosions as nuclear or non-nuclear* 

NRDL has pursued the problem of classifying a reported detonation 
by locating the scone by moans of cm appropriate search based on the fix 
information received from hydroacoustic and/or other coiTolated equip- 
ment, obtaining samples if possible and analyzing pertinent character- 
istics of the area* 

The classification problem splits into two distinct categories re- 
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fleeting two completely different physical characteristics of the scene. 

In the first case, the iinderwater nuclear explosion vents to the surface 
and scatters the attendant radioactive debris on the surface and probably 
into the air. In the second case the radiation is contained below the 
surface due to a combination of depth of explosion and yield of the 
device. This thesis is devoted to analyzing the first of these two 
categories . 

2. Statement of the problem. 

The objective of this study is to determine the optimum searching 
procedure for locating the site of an underwater explosion which has 
vented, in order to classify it as nuclear or non-nuclear. 

3. Means of identification. 

The most positive means of confirming that an explosion is nuclear 
is by a radiochemical analysis of a sample of the affected area. 

After an explosion has been detected and a radiating sample obtained, 
it remains to distinguish the origin of the radiation. In spite of a 
radioactive sample a nuclear detonation could be disclaimed and the 
activity attrlbvrbed to one of the following: 

i 

a. Residue from previous explosions, 

b. Waste dumping, 

d. Reactor discharge, 

d. Accident. 

Radiochemistry can detect fission fragments and normally determine 
the time of detonation by extrapolating back from ratios of certain iso- 
topes present in the sample. This information correlated with seismic, 
hydroacoustic, air acoustic, and/or electromagnetic information of time 
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and location of detonation can readily refute a claim of residue from 
previous explosions or waste dumping as the sovirce* 

Normal reactor discharge does not contain fission fragments and 
would easily be repudiated as a cause for contamination. However, leak- 
age from a reactor primary system could be a source of deceiving activity 
in the sample, but would not explain the detected explosion at the loca- 
tion. 

An explosion of a reactor in water would have attendant debris, but 
with the burden of proof to counter such a claim upon the investigator, 
failure to find debris is unlikely to be satisfactorily conclusive that 
such an explosion did not take place. 

Any surface sample is subject to prior masking by the damping of 
fission fragments to confiise analysis. A determined effort to make 
radiochemical analysis of surface samples inconclusive is almost bound 
to be successful in a water contaminated area if the deceiver is willing 
to invest enough effort. Contaminated air samples from venting of under- 
water bubbles are less vulnerable. 

Ibderwater pulsing of the bubble on its route to the surface leaves 
a trail of discrete thin laminas of radioactive material. To reduce the 
possibility of confusion of the surface sample by a masking attempt a 
further effort should be made to sample some of these lower strata. 

Sampling procedures should recognize biological phenomena which can 
corroborate other evidence. Certain algae have a known affinity for 
Iodine 131 which has a short half life and can be identified as a fission 
fragment. Plankton can also produce evidence of an event. If there are 
land and shallow water in the vicinity, molluscs and tropical plant sam- 
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pies could produce valuable Information, 

Of course, specifically, veil equipped search vehicles could arrive 
at reliable conclusions at the scene with pulse height analysis equip- 
ment, air drop water probes and well schooled personnel* However, it is 
not realistic to Imagine planes particularly outfitted for this searching 
chore* Our forces must be more fleocible for the large variety of mis- 
sions they may be called on to perform* It is envisioned that patrol 
craft could be rigged with the detector described in I3SNRDL-TR-525/1/ or 
that bases from which these planes operate could have these meters avail- 
able* However the more sophisticated equipment is not likely to be as 
well distributed* 

Actual meter readings would be of some assistance in correlating 
with the time of the explosion hy determining decay rates j however, con- 
tinued dilution due to wind and sea movements will make analysis of this 
data much less reliable than water samples* 

4* Detection capability* 

It is worthwhile to consider some of the major factors involved in 
the transmission of sound waves in water* 

The sound reaching a receiver experiences three kinds of degradation t 

a. Divergence - For a point source in a homogeneous medium the 
energy is dispersed in a sphere as it is radiated out so that a receiver 
perceives only an amount of energy proportional to the percent of area 
of the expanding sphere that it can observe* 

b. Attenuation - There is a damping of the energy becaxise the pres- 
sure is not in phase with the condensation duo to viscosity, heat conduc- 
tion and molecular absorption (due to magnesium sulphate)* 

c* Anomalies - There are many anomalies which contribute to the 
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peculieirities of underwater sovuid transmission but the principal ones are 
reflection and refraction. 



As Temperature goes up P decreases; 

As Salinity goes up B increases; 

As Pressure increases P increases but B increases more. 

Therefore the velocity increases as temperature, salinity and pressure 
increase. 

Pressure increases with depth but temperature normally decreases, and 
temperature reductions have an overriding effect on the velocity. In the 
Atlantic a temperature of 2®C is reached at about 4000 feet and remains 
fairly constant from then on. Salinity remains relatively constant and 
velocity of sound increases with depth due to the pressure increase. In 
the Pacific the same conditions obtain at about 2000 feet. Near the 
equator this condition of constancy is reached at a lesser depth. 

Near the surface radical temperature changes are not at all infre- 
quent, Surface disturbances of storms mix the water for a considerable 
depth and sometimes send warm water down while the colder water is 
brought to the surface. This produces temperature inversions which have 
a drastic effect on sound transmission. This will be discussed in more 
detail below. 

As the velocity increases the ray (perpendicular to the wave front) 
curves downward. Conversely, if the velocity decreases the ray curves 
upward. 



The velocity of sound in a fluid 




where B = Bulk Modulus 



P = Density 
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Consider the- following ray diagrams where the left depicts the tem- 
perature chamges with depth and the right shows the ray path (6o is 
original direction of propagation). 

TEMPERATURE RAY PATH 




Fig. 1. Effect of temperature changes with depth on sound ray path. 

Sound chamneib occur when there is a minimum velocity at some depth 
with higher velocities above and below. The region of minimum velocity 
is called the axis. Sound originated on the axis and propagated along 
the axis will not diverge spherically but only cylindrically which moans 
that the losses are cut down. Additionally, there is no major anomaly 
to degrade the ware. 

a. Normal Sound Channel 
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b. Pseudo Sound Channel 
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Refraction due to temperature inversion 



c. At depth of 4000 feet in Atlantic and 2000 feet in Pacific: 




3000* 



RAY PATH 




Fig. 2. Sound channels. 

Sound sources do not have to be in the sound channel in order for 
sound to be propagated along the channel. Rays entering sound channels 
at small enough angles (usually less than 15®) with the horizontal may 
be entrapped and then act as if they originated on the axis. 

These few facts may help the reader to appreciate the following 
comment paraphrased from the preliminary NEL report /?/ 1 

Historically it is apparent that even small explosions under ideal 
' conditions can bo detected at ranges of many thousands of miles' , 

whereas, the most powerful explosions under other conditions may 

1 MILES will bo NAUTICAL MILES (6087 feet) throughout this thesis. 
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go undetected. For example; one pound of TNT exploded on the sound 
axis in the open ocean would be expected to produce a stronger sig- 
nal them a 100 Megaton explosion under other extreme conditions. 
Obviously, then, the hydroacoustic detection of a nuclear burst is 
a complex problem and a generalized statement of the capability of 
such a system is not possible. 

The interfaces between the air and the water or the land and the 
water provide a poor match for transmission of sound waves (high coeffi- 
cient of reflection). As a result surface bursts or undergrotind bursts 
are not well detected through hydroacoustic moans Just as underwater 
bursts are not well detected by seismic devices. For this same reason 
stations are not able to monitor inland seas and locations in lagoons or 
with seamounts between the event and the receiver. 

The velocity of sound in water is about 4800 feet a second and it 
is possible that some stations will receive a signal as long as one hour 
after detonation (less than 3000 miles). It will take at least two more 
hours to correlate station information, evaluate, and call an alert. 

There are many sources of confusion in detecting, Ifydroacoustic 
stations will bo unable to distinguish between chemical explosions and 
atomic explosions. Earthquakes and underwater volcanic disruptions when 
detected will bo a concern and frequently non-discemlble from man made 
explosions. The splash of missiles in deep water, ship propulsion noises 
£md environmental sounds characteristic of the ocean can be distinguished. 
In the Joint Committee Hearings /8/ Dr, H, Brown of Lawrence Radia- 
tion Laboratory estimated that there wore 5000 continental earthquakes 
^ per year that could be suspected of being nuclear explosions of between 
one and five kilotons and that only a few percent could be identified as 
earthquakes by present criteria. About 100 disturbances would be sus- 
pected of being nuclear explosions of five kilotons or more. Since most 
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of the continental seismic false alarms are caused by these earthquakes 
it seems unlikely that the number of hydroacoustic false alarms duo to 
natural causes would bo this order of magnitude* Earthquake disturbances 
are not likely to be transmitted with too much energy through the inter- 
face with the water. Volcanic disruptions below the water may be mis- 
leading, but these are not near as numeroiis as continental esurthquakes . 
The principal cause of false alarm in the hydroacoustic network is moot 
likely to be chemical etxploslons transmitted by the Sound Channel mecha- 
nism* 

5. Alternative courses of action. 

With our objective clear in mind - to determine whether a reported 
explosion was nuclear or non-nuclear - our problem is reduced to search- 
ing the area in question to determine if there is evidence of nuclear 
activity. Only by not obtaining evidence of nuclear activity can we 
assume that the disturbance was non-nuclear (ignoring the possibility of 
non-venting nuclear eocplosions for this study), and even then the conclu- 
sion is not absolute. The certainty with which a conclusion can be 
reached depends on the thoroughness of the search. 

The courses of action that are open revolve about the approach to 
the development of a search plan. 

' a. Should the plan bo based upon a random search approach, consid- 
ering the target to be equally likely in some large area prescribed by 
h 3 rdroacoustic fix information? This would permit parallel sweeps through 
the area in an orderly fashion by several aircraft with no danger of in- 
terfering with each other. It would assure a thorough coverage of the 
area because navigation is simple and reliable. 
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b. Should the plan attempt to optimize seeirch effort with respect 
to time? This would permit target distribution to be any shape that the 
physical situation dictates. The pattern would be formed to search the 
areas of greatest probability continuously with the distributicn of the 
fix being the "a priori" function and an "a posteriori" function formed 
after each search. This approach would be best for the following con- 
siderations! 

( 1 ) Search craft may have to be recalled prior to completion of 
the plan or there may be a cost or utility consideration. 

(2) International ground rules may limit time to search a given 
area or it may be desired to spend as little time as possible in a given 
area due to political or other considerations, 

(3) There is the over present marginal case where just enough 
radioactivity is present at one time in sufficient intensity (because of 
yield, depth of detonation, and/or dispersion) to be detected, but which 
will not be there at some later time. 

A disadvantage to this approach to the search plan is that, unless 
target location distribution is in fact uniform, the navigation problems 
reduce the reliability of conforming to a given pattern, and holidays 
are left in the coverage. The probability of detection may in fact be 
less than a random pattern, although the same coverage effort is invested, 
because of those holidays. Secondly, there is a peculiarity in the 
radioactive pool contour under many conditions which reduces the proba- 
bility of detection at certain times. The character of many pools is as 
follows t 

Upon initial venting a high intensity area is established for a 
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certain distance around "grovind zero*" After some time passes many of 
the fission fragments sink, and fringe area activity is diluted beyond 
recognition* This decreases the contour of a given low isodose line* 

Then the particles that are suspended near the center spread out with 
water and wind movaments and the low isodose contours spread out again, 
even beyond the original area* 

It is impossible to fill in details of "how much?" or "how long?" 
or "what intensity?" because all answers depend on the physical condi- 
tions of venting and water movement* The point of mentioning it here is 
to acknowledge that it could affect the capability of our searching 
effort* 

6* Measure of effectiveness* 

This situation is different from many military operations where the 
accomplishment of the objective is imperative, regardless of the cost* 

The. search would be conducted duidng peacetime and there may be concern 
for the expense involved, particularly if many false alarms are encountered* 

There is always a possibility of search vehicles having to be recalled 
because of more urgent commitments, logistic considerations or forbidding 
weather* As previously mentioned, international agreements may limit 
searching time* 

Finally, there is a matter of using consistent logic in approaching 
the problem* As presented so far, and as the physical situation dictates, 
the problem admits the probablistic nature of the fix information* Again 
ignoring the non*-venting situation, there can never be an absolute assur- 
ance that the explosion was non-nuclear, but the probability of this fact 
approaches unity as the search effort approaches an extremely large value* 
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Therefore, the search will terminate, if a pool is not located, at some 
time before the evidence is absolutely conclusive that no pool exists. 
This will always bo the case and the amount of search effort, or synony- 
mously - time, is inherent in the problem and cannot be Ignored. 

Ifeder these considerations the MEASURE OF EFFECTIVENESS is to maxi- 
mize the probability of detection per time of search effort. 

7. list of variables to be considered. 

a. Height of search. This is discussed in detail in Appendix I. 

The conclusion is to fly as low as operationally feasible. A plot of 
intensity versus height in Appendix I shows the loss of intensity for 
altitudes from 0 to 40>000 feet. 

b. Detector sensitivity. Appendix III is based on a detector which 
is capable of distinguishing 0*2 microroentgens at the meter. In order 
to meet this requirement the meter envisioned will be collimatedto reduce 
the effect of cosmic radiation. 

c. Search spacing. In Appendix II the search theory involved is 
discussed, and it is generally concluded that the spacing depends upon 
the search width (which is the effective visibility in the definite range 
law) and the standard deviation of the location errors. 

d. Search width. The radioactive pool on the surface of the water 
coupled with the width of the detector view determine the search width. 
This is covered in detail in Appendix III. 

e. Target motion. The movement of the radioactive pool, which is 
the object of the search, will depend principally on currents. Appendix 
rv investigates the variability of the water movement and concludes that 
the speed of movement is predictable but the direction is random. 
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f. Fix distribution 



g. Probabilitj of detection. In Appendix III the probability of 
detection is computed for the ideal spacing - given the search width. 
Since the search width is not known in a realistic situation an aver- 
aging procedure must be resorted to. This is developed in paragraph 11 
below. 

h. Search aircraft characteristics. 

i. Navigation errors in locating the center of the fix. 

j. Time late at the scene. 

8. Ifydroacoustic fix distribution. 

The location of an explosion will be estimated by correlation be- 
tween hydroacoustic detection stations. Their estimates will depend on 
physical condtions prevailing, but for this study the errors have been 
classed as having a standard deviation of 2.5 miles or 25 miles. The 
purpose of assuming deviations an order of magnitude apart is to see if 
the search plan is sensitive to such extremes. The distribution is as- 
sumed normal and the probability is represented byr 

f(x,y) = f (r) = J__ EXP(- r^) 

( 20 *^) 

On the beisls of this assumption Appendix II develops the ideal 
search spacing as 0.75 (eO^ where E is the Effective Visibility and (T 
is the standard deviation. 

9. Navigation error. 

Where there are no visual navigational aids available an aircraft 
can expect to arrive at a given location in accordance with some proba- 
bility distribution. According to Bowditch A/ celestial errors in air 
navigation of five to ten miles are considered normal for favorable con- 
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dltlons (pg 676 paragraph 2806), Electronic navigational aids, where 
available, yield more accurate positioning. In this model it is assumed 
that the error due to navigation is normally distributed with a standard 
deviation of five miles. The variance is additive to the hydroacoustic 
fix variance since they are both normal and independent. When the fix 
standard deviation is 2.5 miles the total location standard deviation is 
5.6 miles. When the fix deviation is 25 miles the location deviation is 
* 25.5 miles. 

These figures demonstrate that the navigation error has a substan- 
tial effect on the problem when the fix error is small, but it is rather 
inconsequential when the fix error is largo, 

10, Search patterns. 

Appendix IV discusses curreiits in detail and for this study the set 
is considered equally likely in any direction at a given location and 
time. The drift, on the other hand, can be predicted reasonably well, 
Koopnan /2/ in sections 1,6 and 7,3 discusses a non-uniform distribution 
of targets. When the fix deviation is 2,5 miles the situation falls 
closeOy into the category of a circular dlstribirtion of the target .with 
movement at a known speed in a random direction. The theory is to search 
on the peak of the distribution, then outside, then inside, ,., etc. 

This pattern will give the maximum probability of detection per unit time 
for locating a moving target whose Initial fix is not definite. See 
Appendix V for details of this plan. 

When the location deviation is of the order of 25 miles or more and 
the time late is less than two days the target can bo treated as station- 
ary and the search pattern that applies is described by Koopmen fz/ in 
paragraph A set of expanding squares - squares of uniform cover- 
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age - are superimposed in the search pattern. See Appendix VI for de- 
tails of this plan. 

11. Search width. 

This one parameter remains to be determined. Appendix III discusses 
the geometry for determining the search width and it concltdes that it is 
the width of a uniform pool of minimum intensity minus the detector view. 
Iftifortunately, there Is no way to determine the pool size. There are 
really not enough statistical data available to generally characterize 
pool contours for various depths, yields, sea and wind conditions and as 
a function of time. 

a. Large location error. 

As a matter of Interest the following tables have been computed 
to consider the effect of using a search width for one size pool when 
the pool is actually some other size. These calculations are based on 
the estimated probability of detection given by Koopman /2/. 

P = 1 - (l ) exp -nw 
( s ) s 

n = Number of searches made (one for this table). 

w = Search width. 

s = Spacing of pattern. 

As concluded in Appendix II the spacing is .75 and for this 
problem W = E. The standard deviation considered is 25*5 miles. 
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TABLE 1, Analysis of probability of detection for given search spacinga. 
Search spa cing (S) = A. 11 miles - Computed for a 2 mile diameter pool 



Pool Diameter Search Width 

(Pool of uniform minimum intensity) (W = E) 



1.3 

3.3 



4 

6 

8 

10 

20 

30 

Search spacing (S) = 6.86 

2 



4 

6 

8 

10 

20 

30 



5.3 

7.3 

9.3 

19.3 

29.3 

- Computed for a 4 mile 

1.3 

3.3 

5.3 

7.3 

9.3 

19.3 

29.3 



3.3 

5.3 

7.3 

9.3 

19.3 

29.3 



w 


Probability of 


s 


Detection 


.302 


0.036 


.766 


.18 


1.23 


to 

• 


1.692 


.505 


2,16 


.634 


4.48 


.939 


6.8 


.991 


diameter 


pool 


.1896 


.015 


.481 


.085 


.7725 


.18 


1.062 


.287 


1.352 


.394 


2.813 


.7714 


4.27 


.9262 


iiameter pool 


.1495 


.01 


.3795 


.055 


.609 


.125 


.839 


.205 


1.068 


.289 


2.22 


.652 


3.375 


.8513 



4 
6 
8 

10 
20 
30 

Search spacing (S) = 8«7 - Computed for a 6 mile 

2 * 1.3 
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Sqarch spacing (S ) = 10.2 - Computed for an 8 mile diameter pool 



Pbol Diameter 


Search Width 


W Probability of 

S Detection 


2 


1.3 


.1273 


.0075 


4 


3.3 


.324 


.04 


6 


5.3 


.52 


.095 


8 


7.3 


.715 


.161 


10 


9.3 


.912 


.23 


20 


19.3 


1.893 


.564 


30 


29.3 


2.875 


.78 


Search spacing (S) = 11*5 


- Computed for a 10 mile diameter 


pool 


2 


1.3 


.113 


.008 


4 


3.3 


.287 


.034 


6 


5.3 


.461 


.076 


8 ' ' 


7.3 


.634 


.132 


10 


9.3 


.808 


.191 


20 


19.3 


1.678 


.5 


30 


29.3 


2.55 


.723 


Search spacing (S) = 16«6 


- Computed for a 20 mile diameter 


pool 


2 


1.3 


.0783 


.005 


4 


3.3 


.1987 


.0175 


6 


5.3 


.319 


.04 


8 


7.3 


.44 


.073 


10 


9.3 


.56 


.107 


20 


19.3 


1.163 


.325 


30 


29.3 


1.764 


.525 
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Search spacing (S) = 20«4 - Computed for a 30 mile diameter pool 



Diameter 


Search Width 


W 


Probability of 


• 




s 


Detection 


2 


1.3 


.0637 


.002 


4 


3*3 


.1618 


.01 


6 


5.3 


.26 


.029 


8 


7.3 


to 

lf\ 

• 


.05 


10 


9.3 


.456 


.078 


20 


19.3 


.946 


.243 


30 


29.3 


1.436 


.421 



Figure 3 is a plot of the probability of detection for each of the 
search spacings under considerations against the widths of the pools* A 
plot against search widths would bo the same shape displaced to the left 
0.7 miles. Figure 4 utilizes the same data but uses SPACING vice POOL 
DIAMETER as the abscissa. 

In order to consider this nebulous problem more quantitatively, 
assume that any search width from.O to 30 miles is equally likely and 
investigate the mean probability of detection for each search spacing 
over the span of search widths. 



P = 



P = 



P = 



/ p(w)f (w)dw = / exp(-W)]^ 

S) sjj30 

J. jl-«cp(-w) - H exp(-w) 

30 A I S) S S) 

_1 fw+S exp(-W) + S 6Qcp(-W) (W+1 )] 

30 U S) S) (S Jq 

|30 + S eocp(-20) + S exp(-20) (^1 ) - 2S | 

L S) S) ( S J 



_i 

30 



The Measure of Effectiveness is again called upon to decide which 
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spacing affords the best approach. It is desired to select the spacing 
which yields the maximum probability per unit time of search effort. The 
spacing (S) is inversely proportional to the time of search effort. 

Hence, by maximizing pS the criterion previously established will be 
satisfied . 

TABLE 2. Comparison of spacings using the MEASURE OF EFFECTIVENESS (LARGE 



LOCATION ERROR. 






Associated 
Bool Diameter 


Spacing 

(s) 


Ave. ft'ob. of Detect. 

(TJj _ 


Measure of 
Effectiveness (ISS ) 


2 


4.31 


.71355 


3.07 


4 


6.86 


.561145 


3.855 


6 


8.7 


.47058 


4.10 


8 


10.2 


.409 


4.17 


10 


11.5 


.3633 


4.175 


20 


16.6 


.239 


3.965 


30 


20.4 


.183 


3.73 



From the above table it can be seen that pS is maximized aromid a 



spacing of 11.5 miles. 



By differentiating the expression for pS with respect to S and 

solving for the maximum the following is obtained t 

d(TSS)= d ]s4§f axp(-20)4^ exp(-20)(20+1 )- s£H 
dS dsL 30 S ) 30 S )(S 15 J 

d(lJS ) =0=1+exp(»30)+ 1 S exp(-30)+30 exp(-30)+2 exp(-30)->- S exp (-30 ) 
dS . S ) 15 S ) S S ) S ) 15 S 

15 

The following transcendental results: 

S In (2sf+45§4^) = 30 
( 2Si?-l5S 
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S = 11 satiBfies this exprefslon* 

Substituting S = 11 into th« ccpression for(^) 

^ = 4.18 (p = 0.38) 

As a result of the above - based on the very stringent assumption 
that search widths will vary uniformly from 0 to 30 miles (implying pools 
from 0 to 30.7 miles in diameter) - it is concluded that a search spacing 
of 11 miles maximizes the ratio of probability of detection to time of 
search effort. This corresponds to anticipating an average search width 
of 8.44 miles and a pool diameter of 9.14 miles of uniform minimum 

i 

intensity. 

If the assumption had boon that pool widths varied uniformly from 0 
to 50 miles in diameter the conclusion would have been that a spacing of 
25 miles is ideal. It is apparent that the conclusion is very sensitive 
to this assumption. Statistical data are required to provide a firmer 
base for this conclusion, but until data becomes available the assumption 
as stated is considered reasonable, 
b. Small Location Error, 

When the standard deviation of the location error is only 5.6 
miles the effect of the current is relatively significant. After a delay 
of 24 hours the peak of the distribution would have moved out to a circle 
of 12 miles radius (assuming a drift of 0.5 knots). The search path is 
along the peak at the time of arrival (modifying spiral to an expanding 
square). After the first path the, "a posteriori" probability becomes 
bimodal with a circular peak Inside and one outside. The pattern con- 
tinues to search these maxima. 

It is beyond the scope of this thesis to develop a probability of 
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detection statement for this situation, but it is worth considering the 
tendency of the "a posteriori" distribution. As each peak is searched 
the surface of probability becomes flatter until after one pass of the 
area the distribution approaches a uniform one. As the first search is 
being made, however, the situation is not much different from that de- 
veloped in the case of the larger deviation. In the latter case the 
probability peak was continuously being pushed outward. In the situation 
of the moving target the circxilar peak is being split into two concentric 
rings enclosing the original peak. Each resulting ring is searched, and 
the peak is pushed both inward and outward as the search progresses. 

In order to got some feel for the problem (which is more difficult 
mathematically than physically) the same eocpression for the probability 
of detection is used as a foundation for investigating the search width. 

From a geometry standpoint, after an area of a given peripheral 
size is searched, 0.75 would receive the same intensity of search in the 
current problem as was received in the stationary, large deviation prob- 
lem. The remaining quarter would be only half covered but with twice 
the effort . i 




STATIONARY TARGET, 
LARGE DEVIATION 



MOVING TARGET-,— 
SMALL DEVIATION 



n = 1 




s ) s ) 




Fig. 5 Comparison of search patterns. 
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Ifcder these very rough approximations and the assumptions of the 
previous discussion a table similar to Table 2 ceui be computed* 



Associated 
Pool Diameter 


Spacing 

Isj 


Average Prob* of 
Detection (n) 


Measure of 
Effectiveness (t5S ) 


2 


2.02 


*767 


1.55 


4 


3.22 


*70 


2.25 


6 


4.08 


.654 


2*67 


8 


4.79 


.615 


2.95 


10 


5.41 


*582 


3.15 


20 


7*79 


*472 


3.68 


30 


9*61 


.409 


3.93 


TABIE 3* Comparison of spacings using the MEASURE OF EFFECTIVENESS (SMALL 



LOCATION ERROR). 



This indicates that the largest spacing would maximize the apprcoci- 
mate probability of detection per time of search effort provided the 
actual probability of detection is at least proportional to the above 
estimation and that pool sizes are equally likely in size from 0 to 30 
miles* 

After one pass of the area using the doubly expanding square and a 
spacing of 9*61 miles the target location distribution approaches a UNI- 
FORM one* It is obviously not UNIFORM, but the peaks have been flattened 
and multiplied and the continued movement of the current throughout the 
search builds up the probability of the target being in a path already 
swept* The end result is that the unit mass is best considered to be 
evenly spread over the area* 

As a result of this reasoning a parallel sweep search will maximize 
the probability of detection per time of search effort after the first 
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I 



i 



pass. The expression for probability of detection then becomes: 

P = W f or S>W} P = 1 for S<W, 

S 



Following the same approach as in subparagraph 11 .a. it is desired 
to determine the spacing that meets the established measure of effeetire- 
ness - again under the assumption of equally likely search widths (pool 
diameters) over the range of 0 to 30 miles. 



P = -1 / W dw for S>W 

30 y S 



= J. 

30 



«^0 

vr 

2S 

-^0 



= 3i 
s 



We are interested in pS. Since p cannot exceed 1 there is no sense 
in 'making S smaller than 15 because p would remain constant as S was made 
smaller and pS would therefore decrease. Obviously the maxlmian is reached 
when S = W. Therefore this search should be conducted at a spacing of 15 

miles. 



For the second pass the search area should be a square whose dimen- 
sions are dictated by: 

a. Three times the location standard deviation to provide an almost 
certainty that the pool was included in this circle at the time of the 
fix, plus, 

' b. the product of the drift and the time that would elapse between 
the time of the fix and the time of ccmpletion of sweep of the area on 
this second pass. 

The sweeps should start along the side perpendicular to the most 
frequent direction of the c\irrent at 15 mile spacing. 
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12, Aircraft capabilities. 

It is assioned that an aircraft will be available for the search 
equipped with the proper survey meter and water sample floats capable of 
the following t 

a. Maximum range to scene 1000 miles with minimum of 3 hours on 
station available, 

b. Search speed 180 knots, 

13* Search procedure. 

This problem deals with an undervater burst which vents to the at- 
mosphere, All of the effort of the study has been devoted to an analysis 
of activity at the scene, but the study has presupposed that an air 
sampling barrier would bo established downwind from the blast. Any 
samples obtained from this endeavor would not be subject to masking and 
quite informative about the circumstances of the origin of the blast. 

As previously mentioned, correlation of hydroacoustic information 
may take several hours. It is reasonable to assume that as much as 24 
hours may pass before search craft could arrive at the scene. The appro- 
priate search plan should bo followed depending on the physical situation. 
If a pool or suspected pool is located an attempt should be made to locate 
the area of highest intensity (with duo regard for safety of the crow). 
Intensity readings should be taken and correlated with position and time 
as accurate as possible. Relative positions within the pool are more 
important than exact geographical locations. At the point of highest 
intensity water sampling floats should bo dropped for later recovery by 
surface craft or helicopters. 

The search aircraft should remain in the area as long as safety and 
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operational requirements permit so it can direct other searching activity 
for recovering water sample floats. 

Frequent Intensity readings may provide sufficient information to 
establish decay curves to extrapolate back to the time of detonation. 

Groat accuracy would bo necessary in correlating the movement of the pool, 
height of aircraft, location of aircraft over some relative point of the 
pool, etc, to permit much of an estimate of the decay rate. The best 
procedure would be to keep recording the intensity of the area of highest 
Intensity to provide some semblance of a reference. This would corrobo- 
rate hydroacoustic information. 

Normal search procedures should be followed in marking datum and 
search area, and making reports, 

14* Conclusions, 

The problem was to determine the optimal searching procedure for 
locating the site of an underwater explosion which had vented, in order 
to classify it as nuclear or non-nuclear, 

A hydroacoustic network is anticipated to bo the detecting agent. 

Such a network has no ability to distinguish an underwater nuclear explo- 
sion from a non-nuclear one. Classification must bo accomplished by a 
thorough investigation of the scene for evidence of nuclear activity. If 
no such evidence is foiond a probability statement can be made of the 
nature of the explosion, and upon this information a decision to terminate 
or continue the search can be 'reached, 

A MEASURE OF E^'YECTIVENESS was established to bej 
Maximize the probability of detection per time of search effort. 
Several parameters were investigated to appreciate their effect on 
the problem. 
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a. Height * 

It was determined that the maximum intensity reading at the 
meter occtits on the surface. As altitude is increased the intensity 
decreases up to an altitude above 20,000 feet. After that the increase 
is gradual and almost negligible. Of course at such a height the in- 
terest is only academic since very few waves have sufficient energy to 
survive that thickness of air, 

b. Search Spacing . 

The development for search spacing as presented by Koopman /2/ 
was expanded and adopted for this study* 

c. Search Width . 

The effective search width was determined to depend on the radio- 
active pool diameter which was an unknown quantity. Pool diameters depend 
on yield, depth of detonation, and physical characteristics of the water 
and the exploded device* In order to discuss the problem quantitatively 
an averaging procedure was used, assuming any pool diameter to be equally 
likely between 0 and 30 miles. This is a very brash assumption, and the 
results of the study depend heavily on this weighting function* Tfefortu- 
nately, little statistical data are available which might indicate the 
distribution of this function. This assumed UNIFORM distribution of 
weights provides an unexpected conclusion when the criterion is the MEASURE 
OF EFFECTIVENESS mentioned above. If the weight function favored some 
particular pool diameter it would be reasonable to adopt the most effect- 
ive spacing for that diameter; however if there is no bias toward a par- 
ticular width the conclusion is not so evident. Various spacings were 
concluded to be optimum depending on the standard deviation of the dis- 
tributions considered. 
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d. Location Error . 

There are two categories of location error considered, large 
errors with standard deriatioiw above 25 miles dictated a search plan 
that ignored ciirrent drift. For smaller deviations the drift had to be 
reckoned with in designing the seaurch plan. The location error consid- 
ered the effect of both the initial hydroacoustic fix error and naviga- 
tional errors of aircraft arriving at a given spot in the open ocean. 

When current movement was considered (with small standard deviations 
of location) the direction was assvmed to be random, but the drift can be 
predicted within an acceptable margin. 

e. Search Plans . 

When the standard deviation of the location distribution is 
large the squares of uniform coverage maximize the probability of detec- 
tion per time of search effort. On the other hand if the deviation is 
small, a retiring search maximizes this criterion for the first pass. For 
subsequent passes of the small deviation problem the ”a posteriori" 
rippled distribution approaches a UNIFORM^ and a parallel sweep search is 
considered best. 

Finally, although this comment is beyond the scope of the stated 
problem, it seems worth mentioning in passing that veiT- little special 
equipment is required to perform this function of classifying detected 
disturbances. Water sample floats which can be dropped from an aircraft, 
sample the water and seal themselves - to be recovered by surface craft 
or helicopters, and the survey meters described by NRDL /l/ are the only 
non-standard items needed. These items can be based ashore and on 
carriers and need not bo on aircraft until a plane is assigned a particu- 
lar task of classifying. 
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APPENDIX I 



HEIGHT OF SURVEYING AIRCRAFT 

In order to determine an appropriate search plan the search height 
must be determined. According to reference /l/ the meter will be con- 
structed with a shield forming a cone about the axis to prevent excessive 
interference from cosmic radiation. This shielding forms a central angle 
of 150®. 

As the search craft increases its height the area seen by the meter 
will increase and one would expect this to Increase the intensity at the 
instrument. However, the attenuation due to the atmosphere will decrease 
the intensity at the meter. 

Model assumptions t 

1 . The pool is UNIFORM and INFINITE in diameter, 

2, Gamma tnergy at the surface is 1 MEV eird the '♦build -up" 
factor is not considered, 

3* Absorption is linear varying at a constant rate with dis- 
tance from the source, 

4* Height of meter is much greater than the diameter of the 
meter opening, 

I = Intensity of radiation at the meter, 

A = Area of meter opening, 

I^= Intensity of point sources at the surface, 
h = Height of Instrument, 

u = linear coefficient of absorption in air (2,265 x 10 "^/FT), 

X = Distance of meter from radiating source, 

r = Distance of source along water from point on surface directly 



below aircraft. 



0 = Angle of rotation of ”r” on surface, 

0 = Angle that point radiating source makes with meter, 
w = Solid angle subtended by meter opening at radiating source. 
The portion of the radiation that arrires at the meter from each 



divided by 4 . 

W A cos 0 

4 "" 4 

This assumes that h is much greater than the diameter of A; there- 
fore the projection of the area of the meter opening on a plane tangent 
to a sphere expanded about the point source is approximately eq\ial to 
the projection onto the sphere. 



point sotjrce is equal to the ratio of the solid angle (w) subtended 




Fig, 7 Parameters of integration. 
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r = 



cos 0 = ii 

X 



(1 ) I = 



h s6o 0 



dr = X 



I/^ A h sjE^-xD^dx 



It would be of interest to find a height which would maximize the 

intensity at the meter if such exists. 

The integrand in (1 ) is continuous except at o and the derivative 

can be investigated by applying Leibnitz's Rule* 

( 2 ) ^ = I^A h exp(-u h sec Cf) sec ^ - h eocp(^ ±1 
dh 2 L ^ 

h sec 0 '~j 

ey) (-wx)d^ 



Consider the last Integral in (2), 
h sec 0 






(-tix)dx= 



?h sec 0 

-egro(*-ux)-u ( «c 2 (-ux)dx{ 



"“jh sec 0 

-e 2 C 2 {-ux )-u(ln x«^+ u?x^ - u^x^- t -u4x^ . , 

X 2*21 3*31 4*4! 

,-h 



This term by term integration is justified by virtue of uniform 
convergence of the infinite series expansion, 

I 

Substituting this in (2)i 

( 3 ) dl = -In sec 0 + u h(sec 0-1 ) - u^h^ (sec^ 0-1 ) 
dh 4 

+ u^h^ (sec3 0-1 ) - u^h4 (sec^-1 ) 

18 96 

Sec 0 = 3.8637 
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dl was programmed on the 1604 computer for a large number of values 
dh 

of h. The results showed that the derivative was always negative, be- 
coming smaller in absolute value as height increased. This indicates 
that a plot of the intensity would be an eocponential decay with its maxi- 
mum value at h = 0. 

Assumption 3 in the previous model is not realistic because the 
attenuation varies with density. Since air density varies with height 
the coefficient of absorption must also vary with height. The mass co- 
efficient of absorption is equal to the linear coefficient divided by 
the density and it does not vary with the chemical or physical condition 
of the medium. For air this quantity is 0.0274 cm^/gm. 

According to the Barometer equation derived in STATISTICAL 
MECHANICS I 

P = Pq exp(-mgh) 

KT 

P = Density at height of instrument, 
m = Mass of a molecule of the medium, 
g = Acceleration due to gravity, 
h = Height of the instrument. 

K = Boltzman’s constant. 

T = Temperature (degrees KELVIN). 

P^= Atmospheric density at sea level. 
u*= Mass coefficient of absorption. 

X = Slant range from radiating source to meter. 

0 = Limit of shielding (75®)* 

Returning to equation (l ) and substituting: 
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“ axp[^*RcJ = expjju*xPQ exp(" <agh )J 

KT 

'h sec 0 



I = 



A h exp -u^xP- expC-ggh) jdx 
2 KT J 



( 4 ) 1 ^ = Qxp[ ^*Pi^ exp(-s^h sec ^ - Se(^*Po ®cp(-Sgh)h] 
B dh hsec^P ’■ 



KT 



KT 




h sec 0 



h mg exp(-fflg^)u»Po axpH^^Po exp(-mgh)^ dx 
KT KT ^ X ^ KT ^ 



h sec 0 _ 

+ I gxp(" Pgh) x Jdx 



B = loA 
2 



The last integral of equation ( 4 ) becomes 1 



(5) 



- axp(»ax) 

X 



“|h sec 0 /** h sec 



exp (-ax) dx 

X Where a=u*P^ eocp(-mghJ 



J h h KT 

The first part of ( 5 ) balances the first two terms of (4)* 

/ h sec 0 

exp(»ax)dx 

X 

XI 



( 6 ) 



= a(h mg -1 ) 
KT 



In sec 0 - ah(sec 0-1 ) + a^h^ (sec^ 0-1 )..• 

2*21 J 



Consider the following t 

”a" is positive for all values of h| 

The infinite series is always positive^ 

h mg - 1 is 0 when h = KT 
KT mg 



Therefore there is a maximum or minimum or flex at h = KT 

mg 



d% 

dh^ 






h = CT 

ag 



= b 4 In sec 0-a h(sec 0-1 ) + a^h^ (sec^ 0-1 ). > . 






2 * 2 ! 
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Again "a", B and the series are positive. The conclusion is that 



h = KT is a minimum point, 
mg 

Equation (6) was programmed out to 30 terms for values of height 
from -100 to 40,000 feet and several sensitivity checks were made. 



gated and it was varied from 300®K (81 ®F) to 273®K (32®F). The values 
varied slightly but the structure of the curve remained the same. A few 
points were checked with the temperature at 243 (-40®F) and agedn the 
height of the minimum decreased slightly bui the structure of the curve 
remained the same. This covered the range that average temperatures 
could reach and the minimum varied from 21100 feet to 26070 feet. 

Variations in surface air density (Pq) were looked at but their 
effect were even less significant. 

Variations in ”g" were not looked at, but they are obviously of no 
conceni. 

Lastly, the evaluation was made for 10 terms of the series ard then 
for 30 terms of the series. In both cases the values were the same out 
to eight significant figures. 

As a final investigation the expression for the intensity was pro- 
grammed and curves obtained are presented in figure 8. It is satisfying 
that the same MINIMUM points were reached by two independent approaches 
and two different programs. 

The development was as follows: 



Ae temperature is recognised as a variable over the range invest!- 



/ 



h sec 0 




I = 



h 
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= InAh -exp (-ax) 


h sec 0 i^h sec 0 


-« 1 exp(-ax) dx 


2 L, X 


J ^ - 

h 



= - I/Jthl exp(~h sec 0 a) ~exp (~ha) +a In sec 0-a hfsec 0-1 I 
2 [h sec 0 h J L ^ — J 

= loA I /l -1 V hafln sec 0f(l-1 \ha(sec 0-1 ) 

2 sec 0] ' 2!’^ 

+f - 1 H)h ^a^ ( sec^ 0-1 ). . «-^/ f 1 - 1 a” ( sec^ 0-1 ) , • 

V 2*2r'5J (•^•hl'+ThvTTl'j 

I = 1^741 +ha (-la sec 0+(l-1,^ ha(aec 0-1 ) Q 

In fig. 8 the bracketed quantity is plotted as f(h) - the "Intensity 
factor". 



CONCLUSIONS 

1 • The greatest intensity at the meter will be experienced when the 
height is zero. 

2. The minimum point of intensity is little more them acedemic 
interest since the assumption of a UNIFORM source of diameter at least 
equal to the field of the meter becomes a stringent condition at 25,000 
feet. The pool would have to be 30 miles in diameter. More directly, 
this problem is interested in maximum intensities. 

3. Operational conditions will dictate the height of search with 
the aim of flying as low as feasible on a prolonged mission. Figure 8 
demonstrates the degradation of intensity at the meter (and therefore 
the decrease in search width) for the height of search selected. It is 
not possible to predict the decrease in probability of detection which 
will bo experienced as height increases because the physical conditions 
of the pool of radiation are not known. Since the meter has a threshold 
intensity which must bo reached before a signal can bo recognized, flying 
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at a greater altitude than operationally necessary may degenerate the 
intensity at the meter so that the threshold is not reached - whereas 
it may have been reached at the lower level. 

4. At the altitudes that are operationally feasible the curve will 
give a reliable indication of Intensity factor (f(h)) for any atmospheriftr 
conditions (neglecting the build-up factor effects ) since it is insensi- 
tive to temperature changes below 3000 feet. 
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APPENDIX II 



SEARCH THEORY 



This appendix is based on Koopman's Search and Screening (OEG Report 



tail as it applies to this thesis*.; 

A. Development of probability of detection of a target of lateral 
range X from a straight aircraft trader 



instantaneous probability density of detecting a target. If s remains 



ability of detection during time t. 

If s is a changing density because of relative movement with respect 
to the observer, and is time - dependent, the expression is more compli- 
cated. The instantaneous density s is inherently dependent on range and 
the fact that it vsiries with time is conveyed by the symbol s^. 



No. 56 ) /2/. An attempt is made to present his development in more de- 




1. Consider the case of "continuous looking" where "s" is the 



constant throughout a search p(t) = l-oxp(-st) where p(t) is the prob- 




I C 




( 






Fig. 9 Target's relative track. 



The target is at (x,y) at time t and x = x(t), y = y(t). Let the 
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initial position at t = t* bo = X(t*), - 7(t*), and the final posi- 

tion at t = t" is X^ = X(t"), 7 -^ = 7 (t") with the target tracing a track 
along C in figure 9 above, 
s = sjjx 2 (t)+ 3 T^(t))^J= 

The probabilit 7 of detection along path C isi 

Pc = ^ - “pf-X * 8 ( (x2 ( t )4y^ (t i]‘^)dt^ 



The integral is a line integral along C| if ”w” is the relative 
speed and ” 1 " is the arc length along C from X^^, Y^, then 

Pq = 1 - expX js(Qc^(t)+y^(t^‘^) ^ dt = dl 

^ c V w 

Following Koopman*s notation the exponent is called the sighting 
potential f|c^. This quantit7 possesses the propert7 of €uiditivlt 7 so 
that where C = + C2> then F^C"] = F^^C^ 

In the case where X is the lateral range, X is consteuit and y = wt 
(w is a constant). 





i( |x^-H#^^'^)d1 




^ s(^^4y^^)d7 



According to the inverse cube law of sighting where h << r* 
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Fig. 10 Detection at fixed speed and course. 

When the sweep is "relativelj" infinite 6* = 9" = 90° and - 

( 3 ) M 

uTyi- 

When the observer and the target are on straight courses at constant 
speeds for a long time before and after their closest point of approach, 
the probabilitj of detection is a function of the lateral range X. 

2. Parallel sweeps - 

A target is presumed to be at rest on the ocean in an unknown 
position, all equal areas having the same chance of containing it. A 
search is made along a largo ( ”inf inite " ) number of parallel lines a 
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conmon disteuice apart. This distance is referred to as Sweep Spacing (s) 
eind the problem is to determine the probability of detection of a target, 
P(s), by this configuration. Other questions concerning probability of 
detection, if target speed and course are known or have a distribution 
of their own, can follow in development. 

Consider the lateral ranges of the target from various sweeps 



I 

I 

1 1 1 1 — . . A” 

-Z5 “5 0 5 - - - 

Fig. 11 Parallel sweeps. 

For paths to the left the distances are X, X + s, X + 2s... and for 
those to the right s - X, 2s - X, 3s - X... Both cases can be combined 
in the absolute eocpression 

lateral range = j X - n® | where q ^ X S 



n = 0, ±1, ±2, ±3, ... 

The probability of detection by the lateral sweep when the tar- 
get is at (X, y) is: 

The probability of no detection by the n"^^ sweep is 1 - Pj^; and no 
detection by any sweep is the infinite product IT (1 - Pjj) values 

of n. The probability that at least one sweep detects a target at 
(X, y) is ^ 

( 4 ) <p(x,5}J 



(5) 



([ 1 ^-”^ ‘ 

(X5 

where ^ ~ 
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Consider the lateral distance of the target X as haying a uniform 
distribution between 0 and s. P(s) is the arerage of P(X, s) oyer all 
values of X in the interval* The expected or average value of the 
probability of detection for a given s can bo expressed as a function 
of 8* 

(6) P(S)- i ^ ^ 

This is the solution to the problem but it remains to express it 
more quantitatively* 

The effective visibility (E) is defined as half that sweep spacing 
for which the probability of detection by parallel sweeps is i*e* 
P(2E) = 

In the definite range law detection occurs if smd only if, the tar- 
get happens to be within the definite range R of either of two adjacent 

sweeps* The chance for this is gR = W when s ^2R = W, and unity when 

s s 

s ^ W. From this it follows that E = W (-since W = *5j iL = .5; W = E)* 

8 2E 

For the inverse crube law when h <f< r wo get from equation (3 ) t 
F(X) = 2 h k t and from equation (5)f 0(X, s) = 2 h k 

rx^ r 

This can be restated as 0 (X, s) = 2 k h 2- csc^A* because of the 

W“ s^ s 



following development t 

a. The infinite product expression for sin 3 is 

sin 8 = e TT J_, ) 

b. TAke In of each ~ 

In sin 3 = f 
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c. Differentiate both sides - 

gQS ^ 
sin S 



f -I 

= X f % 

-2 M-'i '- 



d. Which may be rewritten t 

cot 8 = -i - 'S: -4-—X 

Note: This is the Fourier expression for cot 8 and could have been used 

as the starting point of this development. 

e. Break up d. and solve for A and B by partial fractions: 

cot 8 = ( 4?- t V 

From partial fractions = 

A(ntTt 8) + B(nr- *) - 28 
A8 - B8 = 28 
A = 2 + B 
An lT+ BnTT= 0 
2 + B + B = 0 



• t 



c:^ 



(7) Co+ € = -k ' ( 






B = -1 

A = 2 + B = 1 






(8) C ot ^ 'Y ^ JL. 

(9) cot^r ^ 

The only way to see the above steps is to substitute values for n in 

(?) and see that you got two terms for each n which can be represented by 

the expanded limits of the expression in (8). Note also that when n = 0 

in (8) the missing term is 1 which has thoughtfully been provided. 

8 



f. Differentiate (9): 
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Z <r _ ..L — 

-3 — ^ /t.. V i. 



( 10 ) c$c ^ (^s-MirV 



g. Let 3 = 1TX and rearrange (10) 

H- ^ cS 7~^ 



U 



(11) Then C ^ -a. = ^Jl^fcsc^ 



Inserting this in (6): 

1 



= s I" '—'fi 



There is one more bit of wizardry which Koopnan applies 

^ ^ . T5‘ %: J ^ 



Let B = 



) ^ JQ- ~ I 

Tf ^ 



0(B) = P( r ^ ^ J“ *’ (^-^d ojj 



Differentiate with respect to B 



0^B)= eyfC-e <T4C^<S) J& 

= -^ S.yf(-Scoi%) csc'^s ^(5- 

Je>^ 

csc^O dO = -Ki(cot ©) 



when 0 = 0^ cot © = 

© = cot © = 0 



Let cot^© = 



d cot © = dx 



0 '’(B)= f zy.f (- S x^)^K 






■ ^ «K-wt 

Integrate ^ (B) with respect to (B) from 0 to B^ (a particular value 



of B) 
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Let = B 



2 X dx = dB 



X = when B = 0 X = Oj 

, when B = B% X = (b'' )'i’ 

^ 6 r f 2 yJx 

\ CsttT'*- Ory’-Jo y 

(12) 0(B) = ■'^X 

(B^ )^ = X^ = a particular value of X. 



By definition (12) is the error function (erf) based on the normal 
distribution. 

Hence P(S) = 0(B) = erf (b)'^ = erf V'/2kh)'^ 

B\\r' 

For the inverse cube law the effective sweep width W = 2/ 2 kh ^*^ 

\\r / 

(13) P(s) = erf 

s 2 



P(2E) = i = erf W 

4E 

Since erf 0,477 =0,5, then (W)hl = ,477 

4E 

and W = 1 ,076E, hk = ,046E^ 

W“ 

Substituting this in (13) 

P(s) = erf (0,954E ) 

8 

It is of particular interest at this time to note that the proba- 
bility of detection is not particularly sensitive to the law by which 
the spacing (s) is determined - particularly at lower probabilities. 
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12 Probabilities with parallel sweeps. 

Recall the following equations: 

p(X) = 1 - exp^F(X)J for lateral range distribution^ 
p(X) = 1 - expr »2hk 1 for inverse cube law ' 

L w=x2J ^ 

W = 1 .076E from the error function development ' 
hk = 0.046E^ from the error function development 

vr 

Substituting (l5) into (14): 

P(X) = 1 - exp|^.092^^^ 

Which is equation (1 ), the immediate object of the above expansions. 
Before proceeding, let us review the principal assumptions upon 
which the "error function" developient is based. 

a. Target at rest in an unknown position with lateral distribution 
UNIFORM between sweeps. 

b. Observer is on a straight course at a constant speed for a "long" 
time before and after the closest point of approach to the target. (Sweep 
is "relatively infinite"). 

c. Inverse cube law of sighting applies and h\<r« (”h” is the 



Fig. 

(14) 

(15) 



vertical height of the observer and "r" is the horizontal distance from 
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obserrer to center of the target.) 

d. Sweeps are independent. 

B. SPACING (S). 

' Consider the SQUARE OF UNIFORM COVERAGE search pattern (section 
7.3.2 of reference 2). 



Fig. 13 Square of Iftiiform Coverage. 

Let the normal density function ^ 

(16) ^ f (O 1 

represent the probability that the target at the time of the fix be in 
the small region dx dy at the point (x,y), a distance of r from the origin 
and where (T is the standard deviation of fix accuracy. The first problem 
is to find a spacing so that on the initial square ^the probability of 
detection per unit time shall be a maximum. 

As was indicated before, the spacing is not sensitive to the law of 
detection, and the inverse cube law development will be followed so that 
the probability of detection of a target of lateral range x from a 
straight aircraft track is 



^5 S 2S 4S 

CD 
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Reflect for a moment on our problem and the assumptions of the pre- 
vious section. The same error function result is being used. 

1 • The lateral distribution is certainly no longer uniform between 
sweeps; however vdien the standard deviation is large compared to the 
spaoing the approximation of uniformity between sweeps is not very in- 
accurate. 



2. The sweeps are not ’’infinite” in any one given direction, but 
except in the center the length of sweep is relatively leo*ge, and again 
the practical approximation is acceptable. 

3. The Inverse cube law is assumed in the development with h<Or. 
First of all, as already stated, the spacing is not sensitive to the law. 
Caution must be used in proceeding, however, because although the distance 
to the center of the detected configuration can normally be exp>ected to 

be greater than the height of the aircraft, in the case of an extremely 
weak source approaching the minimum detection capability of the meter - 
when directly over the pool - the range to the center would be less than 
the height and the probability of detection would bo overestimated. The 
Instantaneous probability (s) would equal ^ k^^ ^^^vice kh and the error 



would increase as the ratio of height to horizontal distance increased. 

4. Sweeps are still considered independent. 

The probability of the target lying on the strip parallel to the y 
axis between x and x + dx is found by integrating equation (l6) over all 



values of y and obtaining the marginals 

^ I 



( 18 ) 



The peak of the distribution is located at the center. If an 
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indefinite straight flight is made along the y axis and has failed to 
detect the target the distribution is altered and the differential co- 
efficient of (18) no longer represents the lateral density of the tar- 
gets. By Bayes theorem the ”a posteriori” probability considers the ”a 
priori” probability of (18) above and the ”productive” probability of 
not detecting the target on a sweep through the center: 

1- pU) ^ (See equation (17)) 



Bayes Theorem states: 

Let the union of events B^ be the sample space where B^^ 
(i = 1, 2, k) are mutually exclusive and exhaustive. 

Let event C also be defined on the same sample space 






P(Clg^) 

■inwi ti m nwr n-i i“tr wi — i _ 

i r(c\&i)F(&0 

i = / 



In our problem let: 

C = not detecting target 

Bj= target located in strip (x, x+dx) 

The denominator of Bayes theorem can be expressed in the continuous 
case as 



where y represents any strip . 

This integral is just a number, say K, which actually normalizes the 
numerator so that instead of speaking of the probability of the event Bj 
given event C we speak of the density function 



xunouxuii 



K 



- SO that the new density function is proportional to P(c|Bj)f(Bj) and 
this translates to 
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(19) 






In words “ the probability of not detecting the target given that it 
is in a strip (x, x+dx) times the probability that it is in a strip 
(x, x+dx) is proportional to the new density of the probability that the 
target is in a strip (x, x+dx) given that it was not located on the first 
sweep. 

To find the maximum of the new distribution differentiate (19) and 
solve for the x which makes the expression zero. 



T o.o'1'^i-) . , 



2-^r 






X ^ £ 0.6': ~ (E r)' 

So the maximum probability target position is no longer humped at 
zero. Double peaks now occur with a depression at the origin. However 
the distribution is skewed, and the distance between the first and second 
tracks, D, should be slightly greater than .65(E!T)'2' to cover a larger 
volume of probability on the second sweep. 

The probability of not detecting on the second sweep given that the 
target is in a sirip(^x, x+dx) times the probability that it is in a strip 
(x, x+dx) is proportional to the new density of the probability that the 
target is in a strip (x, x+dx) given that it was not located on the second 
sweep. This follows from the same reasoning used in equation (19) 

(- /x) cA r € P - 0 . 2. O'; 3 r 






Integrating this function for all x and solving for a D which would 
minimize the remaining function we would be maximizing the probability 
of detection on the second pass. Koopman presents the solution, arrived 
at by numerical integration, in section 7.3*2 of SEARCH AND SCREENING /2/ 



54 



as 0.75 (E4“)^. This spacing is then used for seecrching about a point of 



fix. 

C. Search Patterns 

\\ 

1 . Square search for a stationary target (Square of Dhiform Cover- 
age"). 

The details of the theory are described very explicitly in 
SEARCH AND SCREENING /2/. This search assumes a circular distribution 
of target position as a result of a probablistic fix. This pattern 
attempts to give a high probability of detection per unit time at the 
beginning of the search when covering the area where the target is roost 
likely to be. 

2. Retiring square search for a moving target (Section 7.3*3 of 
SEARCH AND SCREENING /2/). 

In this development it is assumed that the direction of movement 
is not known but uniformly distributed (i.e., equally likely in any direc- 
tion), while the speed of movement is reasonably well established. 

The theory is to search the peak of the distribution at the time 
of arrival. It has moved outward from the fix a radius of Dt where U is 
the speed of target movement and t is time late. During the search the 
target continues to move and the pattern is to search the peak, then out- 
side, then Inside, then outside, ... etc. This will give the maximum 
probability of detection per unit time of search effort. 

It is recognized that both of these patterns present operational 
problems in implementing. Expanding searches are always frowned upon 
because of the navigational difficulty in making legs precise. This is 
particularly true when spacing is small - less than 5 miles. 
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APPENDIX III 



SEARCH WIDTH 

DSNRDL Technical Report 525 /l/ envisionsan airborne detection in- 
strument capable of sensing radioactive surface activity on the sea when 
the intensity at the meter is 0.2 microroentgens /hr. At a search height 
of 500 feet the intensity of a broad source at the surface would have to 
be 0.74 microroentgens per hour. The meter envisioned would bo colli- 
mated about a central angle of 15O® to reduce cosmic interference. 

A. For determination of the effective search width consider the fol- 
lowing mod el t 

1 . Physical parameters : 

a. Search height - 500 feet. 

b. Meter described above (i.o., shielded around arc of 15O®, 
and capable of detecting 0.2 micror/hr). 

c. Minimum surface intensity present - 0.74 micror/hr. 

d. Meter system time constant about 2 seconds. (The technical 
report anticipates 0.2 sec for the meter and it is presumed that a re- 
corder, alarm and other auxiliary equipment would make system time con- 
stant an order of magnitude greater). 

2. Assumptions: 

a. Surface pool of radiation is circular and of uniform inten- 
sity. 

b. The radiation on the surface ends abruptly at a definite 

boundary. 

c. Definite Range Law of detection obtains (W = E). 

Due to the geometry the meter will see a circular area about 0.61 3 
miles in diameter (See Fig. 14)* 
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Fig* 14 Meter geometr 7 (side view). 

If a circular area of 0.74 micro-r/hr intensity and 0.613 miles 
diameter were present in the search area, the effective visibility (E) 
would be zero since the plane would have to fly directly over the center 
in order to detect it. 

If the pool were larger, the search width could be considered as the 
diameter of the pool minus the meter visibility (0,6l3 miles). 




Fig. 15 Surface geometry. 

Another consideration in determining the effective visibility (E) is 
the meter system response time. In the IBNRDL report /^/ the following 
standeurd formula is provided. 
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(l ) S = So(l ~eacp(^/:^) 

where S = meter indication 
S^= source strength 
t = time over source 
T = syetem time constant 

Neglecting attenixation and build-up, if the meter were over the 
source for ten seconds the meter would indicate over 99^ of the source 
strength. 



For an aircreift traveling at 180 knots the meter would have to be 
over the source for 0,5 miles. This cuts down the effective diameter 
another 0,1 miles for a pool of two miles diameter and a lesser amoiint 
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For this problem a correction of -0,7 miles is made to pool diameter 
to obtain effective visibility. 

Since the size of the pool is an unknovm in the search several sizes 
will be looked at to determine sensitivity. 

For a search of a stationary object using a Square of Iftiiform Cover- 
age pattern assume a fix error (2 standard deviations ) to be 50 miles and 
an error in navigation of 10 miles in locating point of fix. The total 
deviation would be 25*5 miles. For simplicity we will adopt the follow- 
ing expression for the probability of detection! 



(2) p = 1-^1 tnWj exp Koopuan /2/ pg 113 and 114) 

n = # passes over area; W = search width = E; S = spacing. 




Since it is unlikely that data are available for such small intensities 
(0,74 micro r/hr), it is desirable to consider intensities that are more 
likely to be recorded and scale effects from available search data. 

Consider a circular pool of higher intensity, say 0,05 r/hr. If the 
attenuating distance of two configurations were equal, then for each to 
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Induce the same intensitj at the meter^ the product of each area and sur- 
face intensity must equal the same product of the others, 

( 3 ) X, A, ~ ^2. ^ “ Minimum detectable surface 

'^/o ft' intensity when spread 

y /o” ^ uniformly over 

Ig = Ihiiform intensity of new 
configuration 




Fig. 17 Geometry of detector view of high intensities. 

Area of radiation in view of detector is the sum of the segments* 
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(5) Common chord yieldi* 



<1 
5 1 



( 6 ) 



(7) 



( 8 ) 



£ ^ 

The effective search width = 

= 2(h,->-kJ 

Cos I ; 1 ~ 

w = E = 2.( r; ^ -f- 1 (~^ Cos 

From (4 ) t ^ 




Equations (?) and (8) are transcendental and can be solved graph- 



ically* 





From equation (7) 


Prom equation (8) I 


miles 


rad 


ians 


. . radians 


Bool radius 






^ 


02 


ALL 


cr 


:te 


0 


• 0 
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1 
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.033 
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.005 
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.005 
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.081 
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*5 






.01 


.015 


1 
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.01 
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' .01 
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.01 


.097 


4 
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.073 


.01 


.130 


5 


.01 ; 


.066 


.01 


.162 


10 


.01 ^ 


0 


_ ^01 1 


.327 


.5 










1 








* 

i 


2 


.015 


.0745 




\ 


3 


.015 


.062 






4 

if 


.015 


0 






5 

10 




1 


(Cont'd next page) 
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TABLE 








_JFrom equation (?)_ 




From equation (8) 


Pool radius 
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.5 
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1 
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I .035 
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.04 
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0 
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0 
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5, Data for graphical solution to transcend ental equations. 
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Fig. 18 is a plot of these curves. The points of intersection yield 
the solution for 9^ , ©2 various pool sizes. 

Note the limiting pase in which the pool is so large that the arc 
approaches the chord of the detector's circle. In this situation all of 
the £o*ea is described by the segment of the detector's circle. 

k - ^ (©2 “ ® 2 ) 

©2= 0.082 radians. 

The solutions for ©2 range from .0675 radians (3®51 * ) for ^ mile 
riuiius pool to 0.082 radians (4^42*) for an infinite pool. With these 
figures and for the pools under consideration the effective search width 
varies from 

2 (pool radius + .306 mi) for large pools to 
2(0.5 + .306) for the half mile radius pool. 

Therefore the search width is increased by the diameter of the de- 
tector "look" at the surface (approximately) regardless of pool size. 

Remember that a specific example of a pool of O.O5 r/hr with the 
detector at 500 feet altitiide is the case under consideration. 

There is one very unrealistic assumption in the above discussion. 

The intensity at the meter is not just proportional to the area. 

t 

The distance of the meter from the soxirce is of principal concern when 
considering attenuation. Ignoring the build-up factor the inteiTelatlon- 
shlp is as follows: 



Or(><K. ^ 



(9) 



.2 micro rAr = 

X ^rrx^ 
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(The effect of the Barometer equation is considered negligible - see 
Appendix I) 

Equation (9) would be solved numerically for Ag or more meaningfully 
for the l2 required in the pool to register at the meter, 

A rough approximation for the effect of the different distances of 
the source frcm the meter can be obtained by comparing the intensity of a 
particle at the average distance of each configuration. 

The average distance to the particle in the case of the circle of 
minimum intensity is 

I i X r'<drd9 - _ J f (r\h'^y^rcir9(e 

j f r cir ct& 

(For circle) = 0,223 



In the case of the segment the average distance is{ 

X 2 = ,3178 

lot = meter intensity for circle of minimum intensity 
I 2 = meter intensity for segment 
I-|0 = surface intensity for circle 
I 20 = surface intensity for segment 
^ “ Xo A lo C~ for point source, 

Foi"^ HJTy, 1-^0 o.ow\ 

^ 



In order for the meter intensity to be the threshold intensity in 
the case of both patterns, the surface intensity of the segment would 
have to be 1 0,65 the intensity of the circular configuration to compen- 
sate for the different mean distances. 
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In the previous sealing dlsoussion the oharacteristics deduced would 



pertain to a surface intensity of 0*5 r/hr vice 0,05 r/hr* 
The foHoid.ng table will assist in sealing t 
Altitude = 500 ft; = .223; u = 13.786/ml. 



I-jo(rAr) 


A2(mi2) 


limiting 

9? 


h2(nl) 


X2(mi) 


^20/ho 


^20 


W(mi) = 
Pool Size 
plus 
(h2+r2 ) 


.5 


4.375X10-7 


2.18 


.3064 


.3175 


10.6 


5 


.6129 


.05 


4.375X10-6 


4.7 


.3063 


.3175 


10.6 


.5 


.6128 


.005 


4.375X10-5 


10,14 


.3057 


.3167 


10.4 


.05 


.6122 


.0005 


4.375X10-4 


21,94 


.3015 


.314 


9.8 


.005 


.608 


.00005 


4.375X10-3 


47.73 


.2805 


.311 


9.1 


.0005 


.587 


,000005 


4.375X10-2 


107.72 


.1809 


.235 

j 


1.4 




,00005 1 


!.4S7 



If search altitude is 1000 feet the following table can be usedt 
= .446 







Limiting 

®2 


h2(nl) ' 


1 


^ 2 c /^0 


^20 


W(ml)= 
Peol + 


U^ 

o 

• 


1,003X10-3 


10.5 


.6101 


.635 


39 


1.95 


"^722 


.005 


1,003X10-2 


22.5 


.6001 


.628 


34.6 


.173 


1.21 


.0005 


1,003X10-^ 


48.9 


.559 


.599 


19.9 


.00995 


1.17 


.00005 


1.003 


111.2 


.346 


.565 


10.45 


.0005225 


0.96 



TABIE 6. Search widths for various intensities. 
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Fig, 19 Identification of symbols for preivlous tables. 

Interpretation of above tables: 

For scaling purposes use column I20 ^ the surface intensity which 
is known and from which it is desired to scale. 

For example, if the contour of a pool of intensity .05 r/hr is known 
and it is desired to determine the search width, enter column I20 and 
find .05 r/hr. The search width would be the pool size plus .6122 miles 
when searching at 500 feet altitude or plus 1.1 8 miles (interpolating) 
for a search height of 1000 feet. 

Conclusion 

The diameter of the detector view can be added to the pool diameter 
to obtain search width when scaling from surface radiation more than one 
order of magnitude larger than the minimum circular surface intensity. 
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APPENDIX IV 



EFFECT (F OCEAN CURRENTS 



A, General . 

There are many contributing causes of currents, both external and 

Internal to the water. In the deep open ocean wind Is a prominent ln<> 

fluenee* Persistent winds Interact with surface water, Imparting up to 

two percent of the wind velocity to water movement* The direction of 

the movement Is complicated, with the resultant of the "corloHs force" 

and layer friction balancing the actual wind Imparted force* 

P* H* Kuener In MARINE GE0IX)GY /5/, page 29, Is quoted In partt 

Flow In deep water far from the coast driven by a constant wind 
(l*e*, drift currents) at a constand speed, the motion at the sur- 
face Is directed 45® to the right (in the northern hemlsidiere) of 
the wind* Below the surface the direction grad\ially diverges 
farther and farther to the right of the wind and concomitantly the 
velocity decreases* At a depth called the friction depth, which 
varies with the latitude and force of the wind and does not exceed 
200 m, the current la directly opposite to the surface current and 
has fallen to IS of the surface velocity* If bottom friction Is 
left out of account the average movement Is at right angles to the 
direction of the wind* The frictlonfil force exerted by the %dnd 
and the corlolis force at right angles to the movement are then 
equal and exactly opposed In direction* This state Is closely 
attained In nature* Yet there Is some slight friction, and the 
resulting force plus the "corlolis force" Is equal eind opposite to 
the force exerted by the wind 

When the wind just begins to blow the first movement of the water 
will be In the same direction as the wind, but as soon as motion 
of the water sets in the rotation of the earth makes Itself felt* 

In shallow water, where the friction with the bottom Is consider- 
able, the current will be less deflected by the oorloHs force than 
In the Ideal case* In regions like the North Sea and the Baltic, 
where the winds are variable, where the water Is shallow, and where 
the surrounding coasts are at no great distance wind driven currents 
will generally follow the direction of the motive force fairly 
closely* 

Internal characteristics and variants create an Isobarlc graxiieht 
between waters* Salinity, temperature, pressure, level, etc. all have a 
bearing on currentcj however, considering all open ocean currents, none 
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are stronger than the GULF STREAM which has a velocity of up to 3 knots. 
Most open water movements run between 0.3 and 0.7 knots, but It Is very 
difficult to predict the direction of flow at any given location at ai^ 
given time. 

Bowditch /4/ states: 

Current arrows on nautical charts represent average conditions and 
should not be considered reliable predictions of the conditions to 
bo encountered at any given time. 

B. Analysis of Currents in Sample Ocean Areas. 

1. Sot. 

Enclosed are data on current direction and speed. The first 
eight (8) pages (figures 20 - 27) are plots of the frequency of current 
directions (set) for the specified areas for the months indicated. Data 
wore collected by the hydrographic Office from ship reports. The number 
of observations upon which the graphs are based is indicated inside each 
curve. This figure is followed by a single number representing the per- 
cent of observations which reported no current (no set or drift). The 
area represented, broken into 5® X 5® qiiadrangles , is a ten degree strip 
of latitude (25®N to 35®N) across the Pacific from 150®E to 170®W longi- 
tude (about 1-i million square miles). Choice of the area was influenced 
by the large number of observations available over such a large expanse 
of open ocean. The months plotted were chosen arbitrarily to provide a 
throe month period to observe gradual changes (January, February, eind 
March) and other samples during the yeeo* to show overall variation (Juno 
and September). The source of the information is the ATLAS OF SURFACE 
CURRENTS /5/. 

Such raiKiom behavior lends itself to a probablistic formulation 
rather than a deterministic approach. However the metamorphosis of the 
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frequency plots precludes any distribution other than UNIFORM. This 
assumption is the safest one to make when generalising for largo areas 
and over the whole year. 

A bar cheurb of current persistency derived from Hydrographic records 
and compiled by ^br8. Martha Olson of DSNRDL, San Francisco, Calif, is 
enclosed as figure 28. This chart shows the frequency in large ocean 
areas that the most frequent set in smaller 5® X 5® quadrangles occurs 
for three categories of persistency. For example: In 55/^ of the 5® X 5® 
sub areas of the Northwestern Pacific the most frequent direction of the 
set occurs in less than 25/^ of the observations. In 45/5 of the sub areas 
the most frequent direction of the set occurred in 25^ to 50j5 of the ob- 
servations. Table 7 lists the data from which the bar chart was derived. 



2. Drift. 

Figure 29 was compiled by Mrs. Olson from 1961 Ifydrographlc 
Office Information not yet published. It depicts the frequency of drift 
for various mean speeds of the "most frequent current." It shows that 
in over 80^ of the areas the drift of the most frequent current is be- 
tween 0.4 and 0.8 knots. This graph refers only to the North Atlantic. 

Table 8 that follows lists the mean speed and deviation for various 
sectors at various times of the year. In general the deviation is no 

more than 0.2 knots, and the mean is less than the mean of the speed in 

« 

the most frequent direction (because when the set is in other than the 
most frequent direction there is a tendency for the movement to be less 
defined. 

C . Conclusions . , ■ 

1 . SET is an unpredictable variable at any given time. The safest 
and broadest assumption that can be made is that the direction of current 
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at any given time and at any given place ia UNIFORMLY dietributed. 

2. The strip considered in figures 20 through 27 was sampled for 
the month of January to determine the percent of area under the curve 
within the semicircle centered at the prevailing set* Two thirds of the 
total area is included imder this half of the curve, 

4 

3. No matter what the direction of the set is, the drift is 
reasonably predictable. The variance of the drift is small and in view 
of the time late involved in the problem, the solution is not sensitive 
to such small variances. 
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TABLE 7 



A MEASURE OF CURRENT PERSISTENCY IN FIVE SEIECTED AREAS 

PART 1 





REGION 


Numbar of 

Qnadrangla- 

Months 


Ntnabar of 
Observations 


Numbar of Quadrangle-Months 
in Which Most Fraquant Di- 
rection is Obsarv^ the 
Indicated % of Observations 
25% 25^-50^ 50% 


A. 


Northwaatam 
Pacific Ocean 


80 


9,512 


45 


35 0 


B. 


Northaaatam 
Pacific Ocaan 


150 


12,595 


62 


78 10 


c. 


North Atlantic 
Ocaan 


253 


^184,760 


90 


149 14 


D. 


East Indian 
Ocaan 


59 


1,092 


27 


il 1 


E. 


Northvastam 
Indian Ocaan 


96 


10,236 


28 


64 4 




TOTAL 


638 


218,195 


252 


357 29 








PART 2 








REGION 


% of Quadrangla-Mohths in Which Most Fraquant Dirac- 
tion is Observed the Indicated % of Observutions 
25% 25%-50% 50^^ 


A. 


Northwestam 
I^ciflc Ocaan 




56% 




0% 


B. 


Northeastam 
Pacific Ocaan 




A^% 


52i6 


7% 


c. 


North Atlantic 
Ocaan 




35% 


59? 


6% 


D. 


East Indian 
Ocaan 




66% 


29? 


5% 


E. 


Northvastam 
Indian Ocaan 




29% 


67? 


U% 



TOTAL 39? 56? 55? 



^Additional information on th« following paga. 
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TAB IE 7 - Cont*d 



Quadrangle <41ontht Grouping of data from all observationi over the 
yeare taken diiring the indieated calendar month and In the sante quad- 
rangle of approximately 

Source of Data: Measinrement of lengths and widths of spokes of 

current roses in HO IHiblications, Atlas of Surface Currents* 

t 

Courtesy of Mrs. Martha Olson, 

DSNRDL, San Francisco, Calif. 
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APPENDIX V 



SEARCH PLAN FOR SMALL FIX DEVIATION 
Two other errors in location are considered besides the fix error. 

a. Navigational errors of an aircraft in finding the exact 
location to which dispatched are characterized by a normal distribution. 
In this model it is assumed that the error duo to navigation is normally 

V 

distributed with a standard deviation of five miles. The variances of 
the fix error and navigation error are additive with a resulting standard 
deviation of 5*6 miles C|f(2.5)2+(5)2). 

b. Current flow is assumed to be equally likely in direction 
at a given location and time. The drift is dependable enough as pub- 
lished for the area of interest* 

The first pass of the area should conform to the pattern shown in 
figure 31 . The search is designed to investigate the peak of the dis- 
tribution at the time of arrival and follow up with loops outside and 
then inside. This will maximize the probability of detection per time 
of search effort. Koopman /2/ discusses this plan in detail* 
u = drift » 

V = aircraft velocity 
T = time elapsed since initial detection 
The time to go from A to B in Figure 30 is r-^+r 2 
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Fig. 30 Pattern for retiring search. 

The time for the object of the search to move from circle r^ to 

circle T2 = r2-r-| 
u 

To make this equal to search rehiole's time to go from A to 6 , 

T2 - mr-j , where m = v+u » S imi l ar ly r^ = mr2i r4 = nr^j etc.,.. 

v-u 

To prevent overlap and expand or contract about the loop r^ ^ mr^. 

To make this leg be S miles outside the first leg determine r^ suoh that 

r^+r^ = r5-r^-S, from which - mr^+a where a = VS 
Y u V-U 

The distemoe r^ should equal uT. However since in practice the ex- 
panding square is used to approximate a spiral, a mean r^ is selected so 
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that r-j = 0,9ut. To fiirthor allow for the faot that any change in course 
of the pool, at the saiae speed will reduce the outward component, r<f , is 
selected as 0,8ut. 

= mr-^ 

12 = rnly\ +r'| 

13 =ml2 

= mli^+a (outside) 

% = mL^ 

L7 = 

Lg = ml#7-2a (inside) 

L,o = 

- “^0 

L|2 = niL^^+3a (outside) 
etc. 



A sample search is worked out for the following situation t 
Height of aircraft = 500 feet 

Spacing (S) = 9*61 miles 

Drift (u) = 0.5 knots 

Speed of aircraft (t) = 180 knots 
Time late (T) =24 hours 
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= .8(.5)(24) = 9.6 miles 

a = j_E = 180(9.61 ) = 9.65 
T-tt 179.5 



Total mlleetge for this pattern 460 ml. 
Time for pattern - 2 hre 35 min. 

If this plane onl7 hem 3 hrs on station 
available It should spend the last 25 
minutes sweeping the lower left comer. 
Another aircraft will haye to make the 
parallel sweep.second pass which re- 
quires 1 hr 40 min (300 miles ). 



m = 2+a = 160>5 = 1 .0056 
v-u 179.5 

^ = 9.65 

=19.3 

L3 = 19.4 

L4 = 29.2 

I5 = 29.4 
% =39.3 
I7 = 39.5 
Lg = 20.4 
= 20.5 
1*1 0 ” ^ 

1*11 = 

1*12 = 30.3 

1*13 = 

1*14 = 59.7 
1*15 = 60 
Ljg (end) = 60 

The plan leayes a holiday in the lower left comer of the pattern. 

In order to take advantage of the best information the upper right comer 
should be oriented in the direction of the prevailing current so that the 
starting leg is 45® to the left of the most frequent current. 

Note that the center may be adequately covered before outside areas 
are searched - depending on the balance of the time late and standard 
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deviation of location. In caae this happens the pattern should be con- 
tinued with outside expansion only. In the example cited above the cen- 
ter search is just a narrow turn of 270® to right. 

After this search is completed proceed with a parallel sweep as ex- 
plained in paragraph 11b of the main body of this stddy. 
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APPENDIX VI 



SEARCH PLAN FOR LARGE FIX DEVIATION 



For the oase where the fix is not any more certain than within a 
standard deviation of 25 miles or more, the errors in navigation and 
current movement can be ignored. A set of expanding squares - squares 
of uniform coverage - are superimposed in the sq\iare pattern. Koopnan 
fz/ covers the subject in detail. 

The altit\ide of the aircraft should be as low as operationally 
acceptable and for eui example 500 feet will be used. 

Figure 32 is a drawing of the pattern for the following problem: 
Aircraft speed = 180 kts. 

Endurance on station > 540 mllee 
Spacing = 11.0 miles 
Average search width (W=E) = 8.44 miles 
Coverage factor 0 = WL = 4558 square miles 



Maximum search radius = a = (43*7 miles 
Half length of square (Sj^) = 

Number of passes n = a^S = 2 

2Fw 

= 19.5 
S 2 = 33.8 

Leiigth of a side of square = 2^ 

L, = 39 

Lg — 67.6 



-i)] 
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Number of circuits vithin major square (^^) = 



2S 




The second uni form square is tilted 45® to the first. This pattern 
is plotted in figure 32. 

The average probability of detection ist 




If the search aircraft has more endurance than considered above (3 
hours at 180 lets ~ assuming the search area up to 1000 miles from base) 
the plan should be appropriately recomputed. 

For instance if search area were only 500 miles from base 
Station Endurance 1 500 
Coverage Factor 0 = WL = 12660 
Maximum search area radius a = 55*4 
Number of passes n = 3 



= 19.5 
$2 = 33.8 
S3 = 43.6 



Lt = 39 

12 — 67.6 

13 = 87.2 




The average probability in this case would be: 



P = 1-^1+2|^j^2|^= 0.67 
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